Effects of lattice geometry and interaction range on polaron dynamics 
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We study the effects of lattice type on polaron dynamics using a continuous-time quantum Monte- 
' Carlo approach. Holstein and screened Frohlich polarons are simulated on a number of different 

Bravais lattices. The effective mass, isotope coefficients, ground state energy and energy spectra, 
phonon numbers, and density of states are calculated. In addition, the results are compared with 

OA ' weak and strong coupling perturbation theory. For the Holstein polaron, it is found that the 

crossover between weak and strong coupling results becomes sharper as the coordination number 
is increased. In higher dimensions, polarons are much less mobile at strong coupling, with more 

y-^ • phonons contributing to the polaron. The total energy decreases monotonically with coupling. 

Spectral properties of the polaron depend on the lattice type considered, with the dimensionality 
contributing to the shape and the coordination number to the bandwidth. As the range of the 

£Nj 1 electron-phonon interaction is increased, the coordination number becomes less important, with the 

dimensionality taking the leading role. 

. PACS numbers: 71.38.-k 

^ ■ I. INTRODUCTION 

^— > . 

Ctf . 

Interest in the role of electron-phonon (e-ph) interactions and polaron dynamics in contemporary materials has 
recently gone through a large revival. Electron-phonon interactions have been shown to be relevant in the cuprate 
superconductors through isotope substitution experiment^, recent high resolution angle resolved photoemission^, and 
a number of earlier optical 3 *^ and neutron-scattering^ spectroscopies. In the colossal magnetoresistance manganites, 
isotope substitution also shows a significant effect on the physical properties^. It has been suggested that the long- 
range Frohlich e-ph interaction is relevant in cuprates and manganites at any doping because of poor screening 7 , 
although the spatial extent of the interactions is still part of an ongoing discussion^. 

The polaron problem has been actively researched for a long time (for a review see Refs. l9HT(ilh For weak electron- 
phonon coupling A <C 1 in the adiabatic limit huj/Ep <C 1, Migdal theory describes electron dynamics^. With 

■ increasing strength of interaction and increasing phonon frequency u> finite bandwidthi 2 ^ and vertex corrections^ 
become important. The neglect of vertex corrections in the Migdal approximation breaks down entirely at A ~ 1 for 

■ any value of the adiabatic ratio Tiuj/Ep because the bandwidth is narrowed and the Fermi energy is renormalized 
down exponentially. As a result, the effective parameter Xhu/Ep becomes larg In the strong coupling limit A > 1, 
a canonical Lang-Firsov (LF) transformation can be used to determine the properties of the small polaroni 6 . The 
effective mass does not diverge at a critical coupling but instead increases exponentially. No general solution to the 
polaron problem exists for all values of the coupling constant, except in the infinite dimensional limit where dynamical 
mean-field theory (DMFT) can be applied^. It is the enormous differences between weak and strong coupling limits 
and adiabatic and antiadiabatic limits which make the polaron problem in the intermediate A regime difficult to study 

■ analytically and numerically. 

Several methods exist for the numerical simulation of the polaron problem at intermediate A. They include exact 
diagonalizationi^, advanced variational methods 19 (effects of dimensionality are also discussed), conventional discrete- 
' time path integral quantum Monte-Carlo (QMC) algorithms'^! as well as the newer density matrix renormalization 
, group 22 (DMRG), continuous-time QMC 23 i 25 , LF QMC— and diagrammatic QMC 28 i 29 . The methods vary in accu- 
racy and versatility but none can provide all the polaron properties of interest in the entire space of model parameters. 
For instance, exact diagonalizations suffer from the necessary truncation of the phonon Hilbert space, especially at 

■ strong couplings and low phonon frequencies (even then, the total Hilbert space is huge, reducing the number of 
sites and leading to poor momentum resolution), DMRG cannot easily handle long-range interactions, diagrammatic 
QMC and exact diagonalization are inconvenient in calculating the density of states, path-integral QMC slows down 
at small frequencies, and so on. In numerical analysis of polaron models, a complex approach is needed where each 
method is employed to calculate what it does best. 

In this paper we use one of these methods, the continuous- time path-integral quantum Monte-Carlo algorithm 
(CTQMC) to investigate the effects of lattice geometry on the properties of the polaron from weak to strong coupling. 
This algorithm was developed previously by one of us 23 and was based on the analytical integration of phonons 
introduced by Feynmar. 3 — and on an earlier numerical implementation in discrete time by De Raedt and LagendijkSi. 
CTQMC introduced two critical improvements. Firstly, formulation in continuous imaginary time eliminated errors 
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caused by the Trotter slicing and made the method numerically exact for any strength of electron-phonon interaction. 
Secondly, introduction of twisted boundary conditions in imaginary tim o 23 i 31 enabled calculation of polaron effective 
masses, spectra and even densities of states. Although the method is not ideal (it slows down at low phonon frequencies 
because of the condition on the inverse temperature j3 3> (Tiuj) -1 and suffers from a sign problem at small coupling and 
large polaron momenta) it is quite versatile. In particular, it has enabled accurate analysis of models with long-range 
electron-phonon interactions^*^ and a model with anisotropic electron hopping^ 3 .. In addition, it is currently the 
only method that can provide numerically exact polaron densities of states^ 3 * and isotope exponentsSS* 3 ^. We note 
that the variational approach of Bonca et al. has been used to compute the mass isotope effect of the bipolaron, and 
presumably such a method could be used to compute the isotope exponent for the polaron^. 

In this paper, we develop the algorithm further. Because it is formulated in real space it is ideally suited for 
analyzing more complex lattices than the conventional linear, square and cubic Bravais lattices. All that is needed is 
a redefined set of nearest-neighbor hops (or kinks) by which the QMC process moves the polaron path through the 
space. We will analyze and compare several Bravais lattices: linear, square and triangular in d = 2, simple cubic, 
face-center-cubic (FCC), hexagonal and body-center-cubic (BCC) in d = 3, and simple hypercubic in d = 4. To our 
knowledge, this is the first investigation of this kind in the polaron literature. In section ^ we introduce a model 
of electron-phonon interactions and briefly describe the method applied in this work. In section ITTT1 we describe the 
limiting behavior of the polaron problem on different lattices. In section llVl we apply CTQMC to the lattice Holstein 
polaron and calculate dynamical quantities such as the effective mass, polaron dispersion and density of states. The 
effects of changing the length scale of the interaction are discussed in section Finally we discuss the relevance of 
these results in section fvTI 

II. MODEL AND METHOD 

We restrict the model to have only one electron Wannier state |n) and one phonon degree of freedom £ m per lattice 
unit cell. The unit cells are numbered by the indices n or m. The electron hopping is assumed to be isotropic and 
between the nearest neighbors only. The phonon subsystem is a set of independent oscillators with frequency uj and 
mass M . In the real space representation the Hamiltonian reads 

H = H c + iJph + -ffc-ph , (1) 

with 

H c = -tJ2 C l C n , (2) 

<„„') 

m m 

ffc-ph = - ^2 /m(n)c n c n £ m . (4) 

nm 

Here t is the hopping amplitude, (nn') denote pairs of nearest neighbors, and P m = —ihd/dt; m is the ion momentum 
operator. The spin indices are absent because the system contains only one electron. 

The form of electron-phonon interaction is specified via the force function / m (n). The latter is defined as the force 
with which an electron in state |n) interacts with the ion degree of freedom £ m . In this paper, we consider two types 
of the force function: (i) the short-range Holstein interaction / m (n) — n8 nn ^., and (ii) a screened discrete Frohlich 
interaction 

fM = [(m-n) 2 + l] 3 / 2 CXP { R^T) ■ (5) 

In d = 2 (or d = 1), it describes the screened isotropic Coulomb interaction of an electron with linearly polarized 
lattice distortions. The distortions are polarized perpendicular to the plane (chain) of the electron motion, and their 
equilibrium positions are shifted perpendicular to the plane (chain) by one lattice period. In Eqn. Q, k is a constant 
and R sc is the screening radius in units of the lattice vector. Note that the Holstein model is recovered for R sc = 
and the lattice Frohlich model for R sc — > oo. The interaction (0 was introd uced in Ref. l25l to model hole-phonon 
interaction in doped cuprates. Its various limiting cases were studied in Refs. I32l33l3fil37l 
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The focus of this paper is the effects of the lattice geometry on the polaron properties. Each lattice is characterized by 
its coordination number z and the bare electron dispersion E^ = —tN~ x J2 <n n ,> exp [ik • (n — n')]. (For example, 

for the triangular lattice, Ek = —2t[cosk x a + cos(^k x + ^k y )a + cos(ifc x — ^k y )a\.) In particular, we will be 
interested in comparing Bravais lattices with equal z but different dimensionality, d. We will analyze the cases z = 2 
(linear chain), z = 4 (square in d = 2), z = 6 (triangular in d = 2 and simple cubic in d — 3), z — 8 (body-centered 
cubic in d = 3, simple hexagonal in d = 3 (<|| = tj_), and hyper cubic in d = 4) and z — 12 (face-center-cubic in 
d = 3). Note that the minimum bare electron energy is — zt for all these lattices. 

The dimensionless coupling constant of the electron-phonon interaction A is defined as follows: 

m 

This coupling constant is the ratio of the polaron energy in the atomic limit (i.e. at t = 0) to the kinetic energy of the 
free electron zt. It is also convenient to introduce the dimensionless phonon frequency Q = htu/t, the dimensionless 
inverse temperature j3 = /3t, and the dimensionless force / m (n) = k^ 1 f m (ni)- 

The CTQMC method employed here has been described in detail elsewher o 23 ! 25 so here we give a quick overview 
of the algorithm. The initial step is to determine the effective electron (polaron) action that results when the phonon 
degrees of freedom have been integrated out analytically. The action is a functional of the polaron path in imaginary 
time r(r) and is given by the following double integral 

A[r(r)} = J" J'drdr'e-^ (e^-W-A) + e -^/ 2 -|^|)) $ o[r (r), r(r')] 

+ dTTFTm / / drdr'e-^e-^-)($ Ar [r(r),r(r')]-*o[r(r),r(r')]) , (7) 



$Ar[r(r),r(r')] - £ / m [r(r)]/ m+Ar [r(r')] , 



(8) 



where the vector Ar = r(/3) — r(0) is the difference between the end points of the polaron path. These expressions 
are valid when the condition exp (/3o>) S> 1 is satisfied. From this starting point, the polaron is simulated using the 
Metropolis Monte-Carlo method. The electron path is continuous in time with hopping events (or kinks) introduced 
or removed from the path with each Monte-Carlo step. From this ensemble, various physical properties may be 
computed. The ground state polaron energy is 



lim 
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dA\ 1 /v^ Ar 
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where iV s is the number of kinks of type s, and angular brackets denote ensemble averaging. The number of phonons 
is given by: 



N ph = - hm , ( — 
p^oo [3 \dui 



(10) 



where the derivative is taken keeping Xco constant. The polaron band energy spectrum can be computed from: 

£k — £o = — lim s m(cos(k • Ar)) , 



(11) 



where k is the quasi momentum. By expanding this expression in small k, the i-th component of the inverse effective 
mass is obtained as 
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Thus the inverse effective mass is the diffusion coefficient of the polaron path in the limit of the infinitely long "diffusion 
time" (3. Finally, the mass isotope coefficient, a m * — d\nm* /dlnM, is calculated as follows 



= lim — 
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p^oo 2 ((Ar,) 2 ) 
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Here the action derivatives must be taken at constant A. a m * and effective mass averaged over dimensions are related 
to the critical temperature isotope coefficient, a = —d\i\T c /d\nM, of a (bi)polaronic superconductor as 

1 - ^L) , (14) 



A - [ic 

where fi c is the Coulomb pseudo-potentiali^2&. 

Several of the lattices studied here are not bipartite, and observables will quantitatively change with sign of t. In 
this paper t is always taken to be positive (corresponding to s orbitals etc.) Negative t would lead to a sign problem, 
since the probability of the update is proportional to t. One could probably proceed by mapping the negative t 
problem onto a modified positive t problem with long range hopping. We leave such a study for a future publication. 

III. LIMITING BEHAVIOR 

In this section, we discuss limiting behaviors for the polaron problem to demonstrate the physical differences between 
the two limits. The results presented here can also be used to check the accuracy of the QMC algorithm, and are 
compared with the numerical results in the next section. 

A. Weak coupling 

The weak coupling behavior can be computed using a simple second-order perturbation theory. In such a theory, 
the polaron dispersion is given bjiiS: 

(2 ) _ zXuit 2 1 l/ q | 2 

k k E m /^(o)^v^w r (k,q)' (15) 

W(k, q) = £ k _ q + Hlj - £ k , (16) 

m 

(2) 

where N is the number of unit cells. Thus the ground state polaron energy (at k = 0) is €q = — t[z + XT eo (u))] y which 
defines a dimcnsionless coefficient Te ■ A second derivative of Eqn. (|15|1 yields the effective mass 
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_L_ zXcut 2 v |/ q | 2 2W{ 2 - W[W _ i - Ar m . (a) 

m 0l h 2 N ^E m /^(0) ^ 3 (k,q) - m , ' [ ' 



with W[ = dW/dki etc. In calculating the coefficients T m * one should take into account that the expression for 
the bare effective mass changes with lattice type: m^ 1 = 2a 2 t/h 2 for the hypercubic, face-center-cubic, and body- 
center-cubic lattices, but rn^ 1 = 3a 2 t/K 2 for the triangular lattice. The bare masses are not isotropic on the hexagonal 
lattice, with in plane = 3a 2 t/h 2 and out of plane m^ 1 = 2a 2 t/h 2 , so a little care has to be taken in the following. 
The mass isotope coefficient is defined as 

(2 ) = dln(TO*) = u) d ( mpA , 
m ' dln(M) 2{m 0i /m*) duo \ m* J ' 1 ' 

Substituting here Eqn. (|18fl one obtains to leading order in A 

(2) u>9T m *(u}) 

Finally, the number of phonons associated with the polaron in the ground state is 

^ra^« sAr ^- (21) 

The resulting weak-coupling coefficients T can be computed by integrating over q in the Brillouin zone. For the 
Holstein interaction at ui = 1, T are given in Table Q] Note that for all the lattices considered, r m » and r Q , m » are 

identical in all directions, so only one value is presented. 
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TABLE I: Coefficients of the weak coupling behavior for u = 1 and the Holstein interaction computed by numerical integration 
in Eqns. 115I21H . Note that the mass isotope coefficient in d > 1 has a small negative value. 



Lattice 




Te 






Linear 


0.5366(6) 


0.8944(4) 


0.5355(9) 


0.12523(0) 


Square 


0.3610(7) 


1.0161(0) 


0.3610(4) 


-0.00031(3) 


Triangular 


0.3216(8) 


1.1062(2) 


0.3217(2) 


-0.00903(0) 


Cubic 


0.1533(6) 


1.0229(7) 


0.2300(1) 


-0.04022(1) 


BCC 


0.2411(2) 


1.0378(8) 


0.1808(8) 


-0.07974(4) 


Hypercubic (d = 4) 


0.1591(5) 


1.0168(6) 


0.1591(6) 


-0.04533(4) 


Hexagonal 


0.2170(3) 


1.0536(0) 


0.1928(2) 


-0.04461(3) 


FCC 


0.2856(0) 


1.1220(6) 


0.1564(8) 


-0.11021(3) 



B. Strong coupling 



For very strong coupling, the electron becomes very heavy and almost localized, since the act of hopping from one 
site to another by relaxing the local lattice and then distorting the lattice on a neighboring site is very unfavorable. 
The energy of the resulting small polaron is given as E p = —Xzt. As the effective coupling is reduced, other terms 
due to hopping become relevant. The expansion about the atomic limit in the hopping parameter (which is small 
compared to the polaron energy) may be determined through a canonical LF transformation. For nearest neighbor 
isotropic hopping the strong coupling polaron band dispersion is then given bj^: 

e k = -zt\ + e^'^Ek + 0{t/\) , (22) 



^ -, Sm/m(0)/m(b) 

7=1 — ^jm ' (23) 

where b is a nearest- neighbor lattice vector (it does not matter which one). Clearly 7 varies with interaction and 
lattice type. For the Holstein interaction (_R SC = 0), which is purely site local, 7=1. For the lattice Frohlich potential 
(R sc = 00) on a square lattice, 7 = 0.292582, and for the triangular lattice, 7 = 0.197577. Such a polaron is much 
lighter than the Holstein polaron. For the screened Frohlich interaction with R sc = 1, the triangular lattice has 
7 = 0.710852 and the square lattice, 7 = 0.730176. The effective mass is computed as before, using the dispersion 
from Eqn. 

^expf-^), (24) 

and the isotope exponent is 

a m * « ^ • (25) 
The number of phonons is calculated directly from (-ff p h), which happens to have exactly half the magnitude of 

(ifel-ph), so 

« -f = ^ . (26) 

The remainder of this paper will concern simulations of the polaron problem using the CTQMC. The numerical 
results will be discussed with relation to the weak and strong coupling limits. 



IV. THE HOLSTEIN POLARON 



A. Ground state properties from weak to strong coupling 



The main difficulty in studying the polaron problem is the intermediate coupling regime. For A ^ 1, it is necessary 
to rely on numerical results. In this section, we present ground state properties of the polaron. These may be 
determined from CTQMC with no sign problem. 
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FIG. 1: Inverse effective mass for the Holstein polaron as a function of coupling A on (a) linear (b) square (c) triangular (d) 
cubic (e) body-center-cubic (f) hypercubic (d = 4) (g) hexagonal and (h) face-center-cubic lattices, and Q = 1. Also shown are 
the weak and strong coupling asymptotes as calculated in section HTT1 Note that the effective masses of triangular and cubic 
lattices are almost equivalent for A > 1 (both lattices having z = 6). Also the BCC, hexagonal and hypercubic lattices (z = 8) 
have almost identical curves. It can also be seen that the coincidence occurs well before the approach to the strong coupling 
asymptote, indicating that the coordination number is very important to physical properties for most values of the coupling. 
Increased coordination number leads to more mobile polarons at weak coupling and more localized polarons for strong coupling 
with a very fast crossover between the two behaviors. Error bars show the standard error of the points computed using a 
blocking scheme. 
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FIG. 2: Number of phonons associated with the Holstein polaron for increasing coupling with uj = 1. For most couplings, 
the number of phonons is closely related to the exponent of the effective mass. This is expected in the very strong coupling 
limit where mo/m* = exp(— 7JV p h), and the continuation of this behavior to lower A values is of interest. Again, the speed of 
the crossover and the value of A at which the number of phonons reaches the strong coupling saturation value depends on the 
coordination number, not on the dimensionality. 
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FIG. 3: Total energy of the Holstein polaron for increasing coupling A with Q = 1. Note that the crossover from weak to strong 
coupling behavior is very fast for this quantity on the triangular and cubic lattices, but is much slower for the square lattice. 
In general, the speed of the crossover is quicker for larger coordination number. The curve lies below the strong coupling 
asymptote for variational reasons and also lies below the weak coupling result. 
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FIG. 4: Mass isotope coefficient of the Holstein polaron with increasing coupling and uj = 1. At small A, the a m * is negative 
and small, so the effective mass decreases with increased ion mass. Alternatively, the strong coupling behavior shows a clear 
increase of polaron mass with phonon mass. There is an anomalous bump in the isotope coefficient for lattices with z > 6. 
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FIG. 5: Relative differences of physical observables on the cubic and triangular lattices. The total energy is identical to within 
2% on the triangular and cubic lattices, followed by the effective mass, which is within 20% on both lattices. All quantities 
agree well for A > 1. 




FIG. 6: Variation of physical observables with phonon frequency, w on a linear chain with A = 1. The physics remains 
qualitatively similar throughout, with no fast crossover or transition in any property. 
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In Fig. ^ we show the computed inverse effective mass in the regime between weak and strong coupling. Also 
shown on the graphs are the weak and strong coupling results. Panels are as follows: (a) linear (z = 2), (b) square 
(z = 4), (c) triangular (z = 6), (d) simple cubic (z — 6), (e) body-center-cubic {z = 8), (f) hypercubic (d = 4, z = 8), 
(g) hexagonal (z = 8), and (h) face-center-cubic (z = 12), with the same order of increasing coordination number 
followed for all the graphs in this section. The first two panels are the results for the linear and square lattices, 
which were presented first in Ref. The linear lattice has the lowest coordination number of all lattices computed 
here, and as such has the slowest crossover between weak and strong coupling behavior. The inverse mass diverges 
from the weak coupling asymptote at couplings A ~ 1, but the approach to the strong coupling asymptote is very 
slow. The crossover is still slow for the square lattice. In contrast, the change between the limiting behaviors is much 
faster for the triangular lattice, with much larger effective masses in the strong coupling limit. The strong coupling 
result depends only on the number of nearest neighbors z. As such, strong coupling results for lattices with identical 
z have very similar profiles. For example, comparison with the cubic lattice, which also has z — 6, suggests that 
z is the defining quantity for the A > 1 curve, with only the weak coupling behavior dependent on the lattice type 
(we have already seen how the strong coupling limit of the effective mass is exponentially dependent on z in section 
IIIIBfl . For lattices with even larger coordination number, this trend continues. For instance, the strong coupling 
behavior of the polaron on the BCC lattice with z = 8 corresponds well to that of the hypercubic lattice in d = 4 and 
the hexagonal lattice. Increasing z to 12 (FCC lattice) shows an even sharper transition between weak and strong 
coupling behavior. It should be noted that in our calculation the strong coupling asymptote of the effective mass 
depends exponentially on z. We note that in the results from the DMFT with z -> oo presented in Ref. US the 
ratio X/Oj is kept constant as A is changed, so the phonon frequency is continuously increased and the A — > oo results 
correspond to the attractive Hubbard model and are not directly comparable with our results. The computation of 
mo/m* is limited by its magnitude, since small mo/m* means high effective mass and rare hopping, so the ensemble 
has to be sampled for a long time to get good statistics. The curve is truncated at weaker A for lattices with large 
coordination number to take account of this problem. 

In Fig. [21 the number of phonons in the polaron cloud is displayed. Again, a series of different lattices has been 
considered. For most couplings, the number of phonons is closely related to the exponent of the effective mass. For 
small electron-phonon coupling, the polaron cloud is virtually empty. The rate of increase, and the time taken to 
cross over from weak to strong coupling behavior is highly dependent on the lattice type considered. For example, 
on a square lattice, the phonon number is not saturated to the strong coupling limit until A ~ 2, in contrast to the 
z = 6 triangular and cubic lattices, which saturate at A ~ 1.5. For increasing z the saturation value of A decreases, 
with saturation at A ~ 1 for the z = 12 FCC lattice. It is instrumental to consider the ratio N p h/hi(mo/m*) as a 
measure of the departure from strong coupling. 

The total energy is displayed in Fig. |3J The energy decreases monotonically with increased coupling for all lattices. 
The larger the coordination number, the quicker the crossover between weak and strong coupling results. The curve 
lies below the strong coupling limit for variational reasons, and also lies below the weak coupling line. It should be 
noted that, for higher coordination numbers at intermediate coupling, it is necessary to increase the total warmup 
period. This is due to a meta-stable state corresponding to a path with no kinks, with single kinks created and then 
immediately destroyed. The meta-stable state has the properties of the strong coupling (atomic) limit. Once two 
kinks are created (which can take several thousand steps), the system drops out of the metastable state into the lower 
ener gy m inimum, and the correct form for the total energy is found. Such a metastable state has been discussed in 
Ref. l40t To speed up the computation, the path is started with several kinks inserted randomly. In this way, the 
meta-stable state is avoided. 

The mass isotope coefficient a m * is shown in Fig. 0] Note that for Q = 1, most lattices have a small negative 
isotope coefficient for small A. At couplings A ~ 1, the sign of the isotope coefficient flips. For lattices with high 
coordination number, the curve overshoots the strong coupling limit, before converging on the asymptote. 

In order to determine the extent to which the cubic and triangular lattices are similar, Fig. shows the relative 
differences of physical observables on the cubic and triangular lattices. The total energy is identical to within 2% on 
the triangular and cubic lattices, followed by the effective mass, which is within 20% on both lattices. All quantities 
agree well for A > 1. 

Finally, we demonstrate that there is no qualitative change of the physics with phonon frequency. Fig. [fj] shows the 
variation of physical observables as the phonon frequency is reduced. This has been carried out on a linear chain, so 
that the low omega regime can be reached with the QMC solver. The physics remains qualitatively similar throughout, 
with no fast crossover or transition in any property. 
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FIG. 7: Dispersion of the Holstein polaron in the near-adiabatic (Q — 1,A = 1.4) and antiadiabatic (w = 8, A = 8) regimes. 
Also shown is the result from the Lang-Firsov approximation. The band width remains dependent on the coordination number, 
with the polaron spectra of the cubic and triangular lattices having similar width. Alternatively, the shape of the dispersion is 
clearly related to the dimensionality. As expected, the Lang-Firsov dispersion is quite accurate in the antiadiabatic limit and 
fails in the adiabatic limit. 



B. Excited states 



In this section we compute the excitation spectrum and density of states for the polaron gas. The density of states 
could be used to determine thermodynamic and transport properties of the free polaron gas in a similar way to the 
Sommerfeld theory of metals, although some care should be taken in ID, where there are no polarons at half-filling 24 . 
In most materials, it is expected that the density of polarons is sufficiently low to avoid such a breakdown of the 
polaron picture. 

The spectrum of the Holstein polaron in the near-adiabatic (Q = 1,A = 1.4) and antiadiabatic (u> = 8, A = 8) 
regimes is shown in Fig. \7\ for the square, cubic and triangular lattices. Also shown is the result from the LF 
transformation (Eqn. Q22[l). Polarons on triangular and cubic lattices are not very mobile in either regime, with the 
coordination number clearly playing a major role in the band-width. Alternatively, the shape of the dispersion is 
strongly determined by the lattice dimensionality. 

Figure |H1 shows the density of states (DOS), D(e) — X)k^( e ~ £ k)/W for the Holstein polaron in the antiadiabatic 
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FIG. 8: Density of states for the Holstein interaction in the antiadiabatic limit (A = 8 and u> = 8). Again, the effect of 
dimension can be seen in the shape of the DOS, with the logarithmic divergence (van Hove singularity) clearly visible for both 
lattices, while the band width is most affected by the coordination number. 
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FIG. 9: Inverse mass for the Frohlich polaron with u> = 1 (R S c — > oo). Note that the polaron is far more mobile than in 
the Holstein case. The results for the two lattices are very similar, indicating that in the case of long range interactions, the 
dimensionality of the lattice plays a much more important role in the polaron dynamics, with the specific form of the lattice 
unimportant. 

limit (A = 8 and Q — 8). Again, the effect of dimension can be seen in the shape of the DOS, with the logarithmic 
divergence (van Hove singularity) clearly visible for both lattices, while the band width is clearly determined by the 
coordination number (see Ref. |23j for the shape of the DOS of the cubic lattice, which has a completely different 
form) . 



EFFECTS OF SCREENING 



An important question about polaron properties also involves the effects of screening on the electron-phonon 
interaction. Unscreened e-ph interactions make polarons very mobil o 7 i 32 , which leads to strong effects even on the 
qualitative physical properties of the polaron gas. In this section, only the triangular and square lattices are considered. 

In Fig. the inverse effective mass for the Frohlich interaction is plotted. The first point to note is that while the 
inverse mass in the strong coupling limit depended strongly on the coordination number in the Holstein model, this 
is no longer true for the Frohlich polaron. In terms of the way the polaron is formed, the phonon cloud is spread out 
over several lattice sites. This reduces the dependence on lattice type. 

The computation of the energy spectrum is limited by the bandwidth. For very large polaron energies ej. > kgT, 
there is a sign problem introduced by the average of a cosine term (Eqn. (|ll|l ). This means that the spectrum can only 
be computed effectively for small bandwidths. This is not the case for the pure Frohlich polaron, which is significantly 
more mobile than the Holstein polaron. In order to circumvent this problem and obtain some information about 
long range potentials, the screened Frohlich interaction with R sc = 1 is used, which has some properties of both the 
Holstein and Frohlich polarons. 

The band structure and density of states (DOS) of the screened Frohlich polaron are shown in Figs. I1UI and ITT1 
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FIG. 10: Spectrum for a screened Frohlich polaron with R sc = 1, A = 8 and u) = 8. The bandwidth is much larger than in the 
case of the Holstein interaction. The Lang-Firsov result is also shown (light dotted line). It can be seen that the bandwidths 
of the triangular and square lattices have the same order of magnitude as the screening radius is increased. 
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FIG. 11: DOS for a screened Frohlich polaron with R sc = 1, A = 8 and 10 — 8. Again, the lattice dimensionality determines the 
presence of van Hove singularities. The Frohlich polarons on both lattices are expected to have similar transport properties. 



respectively. Even for such a small screening length, the effect on the spectrum is dramatic. Again, the shape is 
given by the lattice type, but now, while the triangular polaron is less mobile than the square polaron, the ratio of 
bandwidths is much closer to unity than for the Holstein polaron. This is a little surprising, since in this case the 
polaron extends over very few lattice spacings but the effects of lattice are essentially annulled. 

Finally to determine how the relative bandwidths become similar in the Frohlich polarons on the triangular and 
square lattices, Fig. ^| shows the evolution of the ratio of the effective bandwidths as the screening length R sc is 
increased in the near adiabatic limit. The ratio of the bandwidths converges to a constant value ~ 1 with increased 
screening length, indicating that the effect of the lattice type is essentially washed out for long range potentials. 



VI. SUMMARY 



We have studied the effects of lattice type on polaron dynamics using a continuous time quantum Monte-Carlo 
technique. The effective mass, isotope coefficient, ground state energy, phonon number, density of states and polaron 
spectrum were calculated for the Holstein and Frohlich polarons. Several lattices were investigated: linear, triangular, 
square, cubic, face-center-cubic, hypercubic, hexagonal and body-center-cubic. The results were compared with 
analytic forms from weak- and strong-coupling perturbation theory. 

The overriding factor for polaron dynamics in the Holstein model is found to be the number of nearest neighbors. 
For instance, the ground state properties for the Holstein polaron on the triangular lattice (z = 6) compare most 
closely with the cubic lattice (z = 6) and not the square lattice (which might naively be expected since square 
and triangular lattices have the same dimensionality). Alternatively, when spectral properties such as the DOS and 
spectrum are considered, the lattice type is strongly responsible for the shape of the spectrum/DOS, but the band 
width is most closely related to the number of nearest neighbors. This is particularly obvious when the Lang-Firsov 
antiadiabatic limit is considered. 
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FIG. 12: Ratio of the effective bandwidth of the screened Frohlich polaron on square and triangular lattices. A = 1.4 and 
ui = 1. Note how the ratio of the bandwidths converges with increasing R sc . This shows that the effects of lattice type are 
essentially washed out for long range potentials. 



When an extended Frohlich polaron is considered, the picture changes significantly. In that case, it is dimensionality 
which most strongly determines the form of the dynamic quantities (such as the effective mass), and not the number 
of nearest neighbors, since the lattice distortion surrounding the polaron is spread over a number of lattice sites. 
The band-widths of the triangular and cubic lattices converge as the screening length is increased, even for moderate 
length scales of the potential of a couple of lattice sites. 

This work was supported by EPSRC (UK) (grant EP/C518365/1). 
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